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ABSTRACT 

NGC 6337 is a member of the rare group of planetary nebulae where a close 
binary nucleus has been identified. The nebula's morphology and emission line 
profiles are both unusual, particularly the latter. We present a thorough mapping 
of spatially resolved, long-slit echelle spectra obtained over the nebula that allows 
a detailed characterization of its complex kinematics. This information, together 
with narrow band imagery is used to produce a three dimensional model of the 
nebula using the code SHAPE. The 3-D model yields a slowly expanding toroid 
with large density fluctuations in its periphery that are observed as cometary 
knots. A system of bipolar expanding caps of low ionization are located outside 
the toroid. In addition, an extended high velocity and tenuous bipolar coUimated 
outflow is found emerging from the core and sharply bending in opposite direc- 
tions, a behavior that cannot be accounted for by pure magnetic launching and 
coUimation unless the source of the outflow is precessing or rotating, as could be 
expected from a close binary nucleus. 

Subject headings: planetary nebulae: individual (NGC6337) - ISM: jets and outflows - 
stars: binaries 
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Introduction 



To explain the complex morphologies and coUimated outflows often observed in 
planetary neb ulae the presence of toroidal magnetic fields and binary n uclei are commonly 



invoked (e.g., 



Garcfa-Segura &: Lopez 



2000 : 



Soker &: Rappaport 



2OOOI ). However, from 



the observational standpoint it has pr oven rather difficu 



evidence for these theories. Recently, 



Garcia-Diaz et al. 



t to de tect the firm supporting 



(I2OO8I ) have discussed the case 



of the planetary nebula NGC 1360 where evidence of a strong stellar magnetic field has 
been detected and magneto-hydrodynamical modeling has successfully reproduced its key 
kinematic and morphological features. In this work we now explore in detail the kinematic 
behavior and morphological structure of NGC 6337, a nebula where convincing evidence 
of the presence of a close binary nucleus exists. NGC 6337 (PN G349. 3-01.1) appears as 
a thick ring with radial filaments and knots, a faint elliptical shell surrounds the ring or 
toroid and a conspicuous pair of condensations bright in [N II] are located at P.A. —43° on 



opposite sides of the ring. T 



le core of this nebula has been identified from time-resolved 



CCD 
days. 



V-band ) 



Hillwig 



phot ometry (iHillwig 



2OO4J ) as a close binary nucleus with a period of 0.173 



(120041 ) assumes a primary mass of 0.6 M© to derive a mass for the companion 



star M2< 0.3 Mq and a binary separation a < 1.26 R©, these characteristics indicate the 
possibility that the binary core underwent a common envelope phase. 



Corradi et al. 



fl2000l ) 



A kinematic analysis of NGC 6337 has previously been made by 
via two long-slits across the nebula with position angles, P.A. = —39° and P.A. = —75°. 
From these data they interpreted the ring-like structure as the waist of a pole-on bipolar 
PN. They also identified the expanding caps, which they labeled features A and B, located 
to the north-west (NW) and south-east (SE) and detected limited regions of high velocity 
that they tentatively identified as corresponding to a point-symmetric bipolar outflow. 
Here we extend that study by using 12 long-slit positions that fully sample the extent of 
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the nebula and allow us to analyze in detail the strange structure of the line profiles and 
combine them with t he morphological inform ation to produce a morpho- kinematic 3-D 



model with SHAPE flSteffen fc Lopez 



2006|) dhttp:// www.astrosen.unam.mx/shape/ [) of 



the main nebula and its complex outflows that resolve the aspect-oriented complexity of 
this object. 

The paper is organized as follows. In §2 we present the observations, in §3 we discuss 
the results, §4 describes the 3-D SHAPE model of NGC 6337 and finally in §5 we summarize 
the conclusions of this study. 



2. Observations and Results 

High-resolution spectroscopic observations and monochromatic images of NGC 6337 
were obtained at the Observatorio Astronomico Nacional at San Pedro Martir, (SPM), 
Mexico, on two observing runs 2006, July 20 and 2007 June 20 - 23 . These observations 



were taken using the Manchester Echelle Spectrometer (MES-SPM) (iMeaburn et al. 



20031) 



on the 2.1 m telescope in a //7.5 configuration. This instrument was equipped with a 
SITE-3 CCD detector with 1024 x 1024 square pixels, each 24 fim on a side. We used 
a 90 A bandwidth filter to isolate the 87th order containing the Ha and [N II] nebular 
emission lines. Two times binning was employed in both the spatial and spectral directions. 
Consequently, 512 increments, each 0'.'624 long gave a projected slit length of 5'.32 on the 
sky. We used a slit of 150 /im wide (=11 km s~^ and 1'.'9) oriented to a P. A. of —43° for 
the majority of the observations and a P.A. = 0° for one pointing across the center of the 
nebula. All spectra and images were acquired using exposure times of 1800 s. The spectra 
were calibrated in wavelength against the spectrum of a Th/Ar arc lamp to an accuracy of 
±1 km s~^ when converted to radial velocity. Deep images of the field were also obtained 
with MES in its imaging configuration in three different filters: [N II] A 6584 A, Ha;-|-[N II], 
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and [O III] A 5007 A with bandwidths of 10 A, 90 A, and 50 A, respectively. 

We reduced the data using standard IRAljl] tasks to correct bias, remove cosmic rays 
and cahbrate the two-dimensional spectra based on the comparison lamp spectra. All 
spectra presented in this paper are corrected to heliocentric velocity (Vhei)- 

Figure [His a mosaic of images of NGC 6337, namely Ha+[N II] , [N II] A 6584 A, 
and [O III] A 5007 A. Since the original images are deep, with high signal to noise ratios, 
here they are displayed at two different digital contrast levels to show the rich structure of 
this nebula. In Figure [21 the slit positions are indicated and labeled on the [N II] image 
of NGC 6337 taken at SPM. We obtained 11 consecutive positions with a P. A. = —43° 
across the nebula and one position at P.A. = 0° to encompass all its main components. The 
bi-dimensional emission line spectra or position - velocity (P - V) arrays from slit positions 
a - f are shown in Figure [3], and from g - 1 in Figure [H 



3. Discussion 

The images in Figure [T] show the bright main ring, its filamentary and knotty structure 
and the bright, low ionization filaments located to the northwest and southeast of the ring 
and that are referred to here as caps. In addition, the Ha;+[N II] image shows a faint 
elliptical outer shell which is also observed in the [O III] A 5007 A image where it displays a 
brightness distribution that resembles a slight "S" shape, sometimes associated with bipolar 
envelopes. A tenuous, high-speed, bipolar outflow that is clearly present in the long-slit 
spectra (see below) cannot be identified from these images and is revealed by the SHAPE 

^IRAF is distributed by the National Optical Astronomy Observatories, which is oper- 
ated by the Association of Universities for Research in Astronomy, Inc. under cooperative 
agreement with the National Science foundation 
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modeling of the nebula to be associated with some of the diffuse material that surrounds 
the caps, most clearly seen in the [N II] A 6584 A image. 

In Figures 3 and 4 a mosaic of three individual bi-dimensional arrays are presented for 
each slit position. The observed Ha and [N II] line profiles are on the left and to the right 
are the [N II] synthetic line profile produced by the SHAPE model. The observed spectra 
are shown here corrected for heliocetric velocity and the offset in declination is with respect 
to the central star. The heliocentric systemic velocity, Vsys, is —70.56 km s~^ calculated 



using slit g, which pass es through the centra" 



the ve locity derived by 



Corradi et al. 



star, this velocity is in good agreement wi th 



fll988h . 



fl2000h and 



Meatheringham. Wood. &: Faulkner 



The spectra show clearly the signature of a thick ring, or more properly a toroid, in the 
form of symmetric bright knots of emission that appear where the slit crosses it. The toroid 
shows line profiles with a wedge-like shape indicating expansion of the section of the toroid 
facing the star, where line splitting is observed. For slit a, the average splitting in the [N II] 
line emission from the top and bottom sides of the toroid amounts to 39.3 km s"^ or 19.7 
km s~^ of expansion velocity along the line of sight. There is some tenuous material close to 
the inner borders of the toroid that share this expansion, as revealed by the velocity ellipses 
in the [N II] line emission from slits e and j that join the emission knots from the toroid. 

The toroid shows a very slow projected radial expansion, in the order of only 1 - 
2 km s~^ on average. This very slow radial expansion together with the nearly perfect 
circular shape of the toroid indicates that it has a very small tilt with respect to the plane 
of the sky (< 10°) and that we are looking at it nearly face on. A crude estimate of the 
deprojected expansion velocity for the ring yields about 11.5 km s~^, assuming a 10° tilt 
with respect to the p lane of the sky. The ang ular outer radius of the ring is 24", adopting a 



distance D = 1.3 pc (IStanghellini et al 



19931 ) to NGC 6337, its linear radius becomes 0.15 
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pc = 4.66 X 10 cm. These values yield a kinematic age for the ring of 1.2 x 10 years. 

The toroid and its inner region contain high excitation gas, revealed by the presence of 
the He II A 6560 line emission. Positions e, h, and i show a bipolar type structure in the line 
emission of this ion. It is however curious that the regions closest to the star, see slits a, f, 
and g, do not show projected He II emission, as if a cylindrical cavity perpendicular to the 
plane of the toroid is indeed fairly void of material close to the core. A cavity produced by 
an isotropic free flowing wind region would be expected to contain the high excitation gas 
at its edge and observable in projection over the star, which is not the case. The presence 
of this cavity is perplexing and to discern its origin will require additional information and 
modehng (out of the scope of the present work). 

The caps are readily recognized in the spectra as the knotty extensions located 
immediately outside the ring emission regio ns, prominent in slit positions e to i. They are 



reminiscent of FLIERS (IBalick et al 



19931 ) though in this case they seem to be possibly 
related to mass ejections associated with the fast bipolar outflows. The cap to the northwest 
is blueshifted and slit f that cuts across the center of this fllament yields a heliocentric 
velocity Vhei= —44.9 km s~^ with respect to Vsys- The corresponding velocities along the 
fllament and just outside of it on both sides (slits d to h) show a decreasing trend in 
velocity from the northeast end of the cap to its southwest tip with values of ^hei= —47.0 
km s~^ (slit d), —45.1 km s~^ (slit e), —40.6 km s~^ (slit g), and —31.0 km s~^ (slit h) with 
respect to Vsys- This effect might indicate that the fllament that forms the northwest cap 
(and associated surrounding material) is twisted or tilted with respect to the sightline. The 
cap on the southeast is redshifted with respect to the systemic velocity and slit g yields 
Vhei= +64.7 km s~^ for this region. The slits on either side of slit g, slits f and h, provide 
very similar values, Vhei= +64.9 km s^^ and +64.8 km s~^ respectively with respect to 
Vsys- Considering the inclination angles with respect to the line of sight derived from the 
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SHAPE model for the caps, their kinematic ages range from 1 — 2 x 10^ years, i.e. these 
are relatively young structures compared to the ring. The diffuse envelope material that 
surrounds the caps is apparent in the Ha line profiles of Figures 3 and 4 where it is seen to 
follow the same velocity trend as the caps. 

An outstanding feature of the [N II] emission line profiles in NGC 6337 are the very 
high velocity components that run from Vhei= —210 km s"^ to \4ici= + 215 km s~^ with 
respect to the systemic velocity, Vgys, corresponding to bipolar, coUimated, jet-like outflows. 
Slits b to g show the line profiles with the redshifted components (east, southeast side of 
toroid) spanning a range of —275 km s^^ < \4ici < ~ 214 km s~^ whereas slits g to 1 show 
those blueshifted (west, northwest side of toroid) within a range of +40 km s"^ < ^hci < 
+ 145 km s~^; the line profile from slit g shows both, i.e. slit g, that passes close to the 
central star shows both sides of the high-speed bipolar outflow, covering a velocity range of 
over 400 km s^^ . The high-velocity components in all slits display their maximum velocity 
close to the star and this decreases with distance from it, contrary to what is usually 



observed in the so-caf 
from the source (e.g.. 



ed Hubble flows wher e the outflow increases its velocity with distance 



Meaburn et al. 



20081) . 



4. SHAPE modeling 



In order to disentangle th e 3-D geometry and ki nematic structure of NGC 6337 



we used the program SHAPE (ISteffen fc Lopez 



20061 ). SHAPE is a morpho-kinematic 



modeling tool that allows the user to reconstruct the 3D structure and observed spectral 
line profiles using expanding geometrical forms. SHAPE uses as reference monochromatic 
images and observed position - velocity diagrams to reproduce the 3-D structure and 
kinematics of the object. Particles are distributed over a surface, or throughout a specified 
volume, and are assigned a specific velocity law and relative brightness. Several particle 
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systems can be used and each can be assigned different velocity laws to form a complex 
object. The resultant 2-D image and spectral information are then rendered from the 3-D 
model and compared with the real data. 

To model NGC 6337, we used the [N II] image and position-velocity spectra since all 
the main kinematic components are present here. Our model was built with a torus and 
conic surfaces for its knotty structure. Segments of spheres were used to model the caps. 
The bipolar coUimated outflows were built from elongated cylinders, flattened and bent. 
The observations indicate that the projected width of the redshifted section of the bipolar 
outflow is slightly wider than its blue counterpart and it has been modeled accordingly. 
This difference in widths may result from projection effects along the line of sight such as 
a slightly twisted bipolar outflow, as may be indicated by the kinematic analysis of the 
northwest cap (see section 3) however, since a projected twisting effect of the outflows 
cannot be disentangled from the emission line profiles of the coUimated outflows, this 
potential twisting has not been considered in the model. Figure [5] shows the resultant 3-D 
mesh of the model before rendering. 

The results of the final rendered model are shown in Figure [6] where they are compared 
to the integrated observed data. Panel a shows the [N II] SPM image of NGC 6337 and 
panel b shows the sum of all of the observed spectral line profiles for slits b - 1, simulating a 
slit equivalent to the width of the region sampled by slits b - 1, but preserving the spectral 
resolution. Panels c and d show the corresponding rendered image and integrated synthetic 
P - V array from the same slits b - 1. The main structural components are labeled in the 
model image and integrated line profile. The individual rendered synthetic [N II] P - V 
arrays from the model are shown in Figures 3 and 4 next to the observed ones, where it 
is apparent that they provide a reasonable match in all instances. Finally, we present the 
3-D rendered representation of NGC 6337 in Figure [7] displayed at various viewing angles. 
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The first frame is sliown witli tlie nortli rotated 43° counter-clockwise, equivalent to liaving 
slits b - 1 aligned vertically. The next two panels are rotated on the y-axis by 45° and 
90°, clockwise, respectively. As an additional test of the overall goodness of our model we 
have produced synthetic line profile s from slits located at P.A. = — 39°and P.A. = —75°, 



Corradi et al. 



(120001 ). obtaining also a good match 



corresponding to those observed by 
with their observations. 

The present SHAPE model is able to reproduce the basic 2-D morphology and the 
set of emission line profiles that provide a representation of the third dimension of the 
nebula through the radial velocity component. Slightly different geometric forms could 
have been used to build up the final model, but in the end there is only a very limited 
set of solutions that are able to replicate the complex P - V diagrams of this nebula. A 
key advantage of SHAPE for objects like NGC 6337 is its ability to model independent 
structures, each with its own velocity law, and then merge them into a single product. We 
have used velocity laws of the type v = k ■ r /tq where k is a constant, r is the distance from 
the source and ro is the distance at which the velocity k is reached. The values for k and 
ro have been chosen to match the observed velocities for the various components of NGC 
6337 and to provide reasonable distance scales along the line of sight. However, the model 
cannot place restrictions on, for example, the length of the bipolar collimated outflows or its 
detailed structure nor on the thickness of the toroid or the precise distance of the caps from 
the toroid. Nevertheless, the present model allows a good understanding of the complex 
structure and outflows of this object that otherwise are rather difficult to visualize. 



5. Conclusion 



A thorough analysis of the kinematic structure and morphology of the planetary nebula 
NGC 6337 has been carried out. The nebula is composed of a conspicuous thick ring or 
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torus and fast (Vhei ^ 200 km s~^ ) bipolar, collimated, outflows that bend in opposite 
directions in a point-symmetric way. The ring does not seem to be the coUimating agent 
for the bipolar outflows. The torus is slowly expanding radially, at a rate of only a few 
kilometers per second, and also shows an internal expansion of ~ 20 km s~^. It should 
be noted that thick rings or toroids are uncommon as main morphologies of planetary 
nebulae, although equatorial density enhancements are common in axi-symmetric nebulae 
such as ellipsoidal and bipolar. The faint outer envelope of the nebula seems to follow 
the emergence of the poloidal outflows silhouetted in projection and it is likely that this 
envelope material is what remains of incipient bipolar lobes blown at an early stage of 
development of the nebula. In this case only traces of these lobes, close to the ring, remain 
now since no additional indication of their presence is apparent in the narrow band images 
nor in the emission line spectra. The caps may be mass concentrations of this material 
that have been pushed aside by the collimated outflows in their way out or a late bipolar 
episodic event, given their shorter age. 



Hill wig 



(120041 ) for the binary core. 



Considering the current parameters derived by 
it is likely that this underwent a common envelope episode that may have influenced the 
formation of the equatorial density enhancement into a thick ring. The region close to 
the core appears as a cylindrical cavity void of material; it is unclear how this cavity 
may have formed. The rich, knotty structure apparent in the thick ring has most likely 
resulted as a consequence of instabilities produced by the erosive interaction with the 
radiation field and winds from the core t hroughout the evolution of the system. We confirm 



the suggestion by iCorradi et al. 



(120001 ) of the presence of a point-symmetric, bipolar, 
collimated outflow; SHAPE modeling of our data has revealed the structure of such a jet 
for the first time. Although the binary interaction may have spun up the white dwarf 
central star, favoring a magnetic launching and collimation of the outflows, it is unlikely 
that a magneto-hydrodynamic mechanism can be solely responsible of the extreme bending 



Garcia-Segura &: Lopez 



fl2000f ) have 



of the bipolar, point symmetric, colhmated outflows, 
shown in these cases that although point-symmetric structures can be obtained from MHD 
models, the collimated outflows interact with the inner walls of the wind blown lobes to 
form point-symmetric structures, as in the case of the planetary nebula Hb 5. The extreme 
characteristics of the bipolar, collimated outflows in NGC 6337 require them to emerge 
from a rotating or precessing source, provided here in a natural way by the close binary 
core, in addition to a MHD coUimating agent. With all these characteristics, NGC 6337 
can be considered an archetypical example of the potential influence of a close binary core 
on the evolution of a planetary nebula. 
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Fig. 1.— Images of NGC 6337 obtained with MES - SPM in the hght of Ha+[N II] (left 
column), [N II] A 6584 A (center) and [O III] A 5007 A (right), these are displayed at two 
different dynamic ranges to show the rich structure of the nebula 




Fig. 2. — The location of each sUt position is indicated and labeled on the [N II] image of 
NGC 6337 




Fig. 3. — Mosaic of bi-dimensional position - velocity (P - V) arrays. For each slit position 
there is a group of three P - V arrays, the observed Ha and [N II] A 6583 A hne profiles and 
the corresponding synthetic [N II] 6583 A P - V array modeled with SHAPE. Slit positions 
a - f are shown here. 




Fig. 4. — As in Figure 3, but for slit positions g - 1. 
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north 





Fig. 5. — SHAPE mesh model of NGC 6337 before rendering, shown at two different orien- 
tations 
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Fig. 6. — Top panels, frames a and b, show the [N II] image of NGC 6337 and the integrated, 
observed [N II] P - V arrray from shts b - L Lower panels, frames c and d, show the synthetic 
image and integrated synthetic [N II] P - V array from slits b - 1, modeled with SHAPE. 
The main components are labeled in both the image and the integrated line profile. 
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Fig. 7. — 3-D rendered representation of NGC 6337 displayed at various viewing angles. The 
top frame is shown with north rotated 43° counter-clockwise, equivalent to having slits b - 1 
aligned vertically. The next two panels are rotated on the y-axis (defined along the vertical 
direction) by 45° (center) and 90°(bottom) clockwise, respectively. 




